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Selected results obtained using THERMUS, a newly-developed statistical-thermal model analysis 
package, are presented. 



I. INTRODUCTION 
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The statistical-thermal model has proved extremely successful in applications to relativistic collisions of both 
heavy ions and elementary particles (cf. pj for recent surveys and references therein). In light of this success, 
THERMUS 0, a thermal model analysis package, has been developed for incorporation into the object-oriented 
ROOT framework This follows closely on the release of SHARE, a statistical hadronisation package, written in 
FORTRAN and Mathematica Q. 

At present THERMUS treats the system quantum numbers B, S and Q within three distinct formalisms: 
a grand-canonical ensemble, in which the quantum numbers are conserved on average; a strangeness-canonical 
ensemble, in which strangeness is exactly conserved, while baryon content and charge are treated grand-canonically; 
and, finally, a fully canonical ensemble in which B, S and Q are each exactly conserved. A detailed theoretical 
overview of the statistical-thermal model as applicable to THERMUS, as well as a description of the structure of 
THERMUS, can be found in @. 



II. STATISTICAL-THERMAL MODEL PREDICTIONS USING THERMUS 

In a chemical analysis within the grand-canonical ensemble, the statistical-thermal model requires 6 parameters as 
input: the chemical freeze-out temperature T, baryon-, strangeness- and charge chemical potentials fiB, an d fJ-Q 
respectively, strangeness saturation factor 75 and fireball volume V. When restricting one's attention to particle or 
density ratios, the volume dependence largely falls away, unlike in canonical treatments where the particle densities 
within these ensembles differ from those in the grand-canonical ensemble by correction factors which are strongly 
volume-dependent as seen in Fig.^ 
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FIG. 1: The canonical correction factor of a S = +2 hadron in the strangeness-canonical ensemble as a function of the fireball 
temperature and radius, assuming full strangeness saturation and S 'fireball = = — 0. 
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FIG. 2: Contour plot of the final K + /ir + ratio in the T — \ib plane with /is constrained to ensure strangeness neutrality, /j,q = 
and 7s = 1. Also shown as the dashed line is the freeze-out condition E/N = 1 GeV. 
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FIG. 3: The statistical-thermal model prediction for the A/ (n) ratio generated by THERMUS. 



Assuming full strangeness saturation (i.e. 7s = 1) and constraining the strangeness density and the ratio B/2Q in 
the model leaves just T and \ib as free parameters. This allows all density and particle ratios to be plotted in the 
T — i-iB plane in which the freeze-out condition 6], E/N = 1 GeV, defines a curve T(/zb) (Fig. [2J ■ 

From fits to experimental data it is possible to determine the energy dependence of the thermal parameters at 
chemical freeze-out, as done in 0. In this way, the statistical-thermal model is able to make predictions as a function 
of collision energy. As an example, in Fig.[3]the THERMUS prediction for the A/ (ir) ratio, with (tt) = 3/2(n + + 7r~), 
is shown. In another application, one sees in Fig. 0|a breakdown in the feed-down contribution to the 7r + . At low 
energy the majority of feeding comes from the baryon sector, while mesons dominate at high energy. This highlights 
the need for a thorough treatment of resonances within statistical-thermal models. The THERMUS distribution 
includes a text file listing all hadrons up to mass 2.6 GeV with u, d and s quarks listed by the Particle Data Group 0, 
as well as text files listing the particle decay channels. 
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FIG. 4: The decay contribution to the n + from various sources (weak decays excluded). 
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FIG. 5: The centrality dependence of the strangeness saturation factor, 7s, as extracted from Au+Au data at ^/snn = 130 
GeV 0. 



III. THERMAL FITS WITHIN THERMUS 

Fits to experimental data are easily accomplished within THERMUS. Since feed-down corrections and reconstruc- 
tion efficiencies are often particle- and experiment specific, THERMUS is capable of ascribing to each particle yield or 
ratio of interest its own particle set and associated decay chain (or separate sets for the numerator and denominator 
in the case of ratios). Use is made of the ROOT TMinuit fitting class. 

In a recent application of THERMUS ||, particle ratios measured in Au+Au collisions at ^/snn — 130 GeV by 
all four RHIC experiments, and converted to a common centrality binning in |lfj|. were analysed. A grand-canonical 
treatment was employed with quantum statistics and resonance width included for all particles, fiQ = 0, and T, fiB, 
Us an d 7s as fit parameters. While the temperature and chemical potentials were found to be largely centrality 
independent, the strangeness saturation factor, 75, was observed to increase with participant number towards unity 
in central collisions (Fig. |3J. This confirmed previous findings at SPS and RHIC energies [Tol ITU that deviations 
from equilibrium conditions diminish with increasing collision centrality. 
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IV. CONCLUSION 

In conclusion, an analysis package, THERMUS, has been developed allowing thermal model analyses to be 
performed within the ROOT framework. Within THERMUS statistical-thermal model calculations are possible 
within three distinct formalisms with a complete treatment of hadronic resonances up to a mass of 2.6 GeV. Fits to 
experimental data are furthermore easily achieved. 

THERMUS is available freely from http://hep.phy.uct.ac.za/THERMUS/, along with a detailed user guide and 
installation instructions. 
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